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Abstract. We present a photometric and spectroscopic study of the blue straggler S 1082 in the open cluster M 67. 
Our observations confirm the previously reported 1.07 day eclipse light curve and the absence of large radial- 
velocity variations of the narrow-lined star. However, we find two more spectral components which do vary on 
the 1.07 day period. We conclude that the system is triple. We fit the light and radial-velocity curves and find 
that the total mass of the inner binary is more than twice the turnoff mass and that the outer companion to the 
binary is a blue straggler on its own account. We briefly discuss formation scenarios for this multiple system. 
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1. Introduction 

Blue stragglers are stars that are bluer and more luminous 
than the main-sequence turnoff of the population to which 
they belong. Hence they appear to be either younger than 
the other stars or have main-sequence lifetimes that exceed 
those of stars of similar mass. Blue stragglers are easiest 
to distinguish in globular and open clusters but are also 
found among field halo stars. Formation scenarios that 
require anomalous evolution of a single star are considered 
less likely than explanations that involve the interaction 
of stars in a binary or in a multiple-star encounter. For 
example, a blue straggler may be formed when one star 
in a binary accretes a substantial amount of mass from 
its companion or when two stars merge either in a binary 
or in a collision. In the last case, blue stragglers carry 
information about the high-stellar density environments 
in which they are found. Reviews on blue stragglers are 
given by e.g. Stryker (1993) and Bailyn (1995). 

The blue stragglers in M 67 have received frequent at- 
tention. Ten of them were monitored for radial-velocity 
variations for nearly twenty years; this revealed one short- 
period (4.2 day) eccentric binary and five circular and ec- 
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centric binaries with periods between 850 and 5000 days 
(Milone & Latham 1992, Latham & Milone 1996). Leonard 
(1996) and Hurley et al. (2001) concluded that not one of 
the proposed blue-straggler formation mechanisms alone 
can explain the properties of these binaries. 

S 1082 is one of the remaining four blue stragglers with- 
out an orbital solution. Its spectrum shows the lines of an 
early F-type subgiant. These lines show only small radial- 
velocity variations (v ra d <7 km s _1 peak-to-peak, Mathieu 
et al. 1986) and little rotational broadening (i> rot sini =4- 
11 km s _1 , van den Berg et al. 1999). This appears to 
be in contradiction to the eclipse light curve with a pe- 
riod of 1.07 days found by Goranskij et al. (1992). The 
relative velocity v of two stars in a 1.07 day binary is 
w=208 (M/M0) 1 / 3 km s^ 1 where M is the total mass of 
the binary; the stellar rotation in such a short-period bi- 
nary is expected to be synchronised with the orbit, giving 
^rot = (R/a)v, with a the distance between the binary 
stars and R the stellar radius. The eclipse light curve ex- 
cludes the low inclination as well as the extreme mass ratio 
required to bring the observed velocity limits in agreement 
with the predicted values for an F star in the binary. We 
therefore tentatively conclude that S 1082 is a (visual or 
physical) triple. 

Evidence for a short orbital period is also furnished 
by the variation detected in a broad, shallow component 
(Mathys 1991) in e.g. the Ha line on time scales of hours 
(van den Berg et al. 1999). Evidence for a wide orbital pe- 
riod of about 1000 days was found from the radial- velocity 
variations of the narrow lines of the F-star (Milone 1991). 
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Table 1. Log of the photometric observations. For each 
run we give the date of the observations, the telescope and 
filters used and the typical exposure time for each filter. 



Run Dates 



Feb 8-19 1999 

Dec 25, 26 1999 
Feb 13-16, 20 2000 



Telescope 

0.91m ESO 
Dutch 

lm ING JKT 
lm ING JKT 



Filters t cxp (s) 

U B 300 120 

V Gunn i 120 120 
BV 75 30 

UBVI 35030 
15 8 



We have collected multiband photometry of S 1082 to- 
gether with high-resolution spectra that sample the light 
curves at different phases. Our goal was to investigate the 
eclipses found by Goranskij et al., to monitor the vari- 
ability of the second spectral component as function of 
photometric phase and to model the light curve together 
with the radial-velocity curves. Sect. || describes the ob- 
servations and data reduction. The results of the spectro- 
scopic and light curve analysis are presented in Sect. || and 
Sect. |. In Sect. | we summarise our interpretation of the 
nature of this blue straggler. 

2. Observations and data reduction 

2.1. Photometry 

S 1082 was monitored in the U, B, V, I and Gunn i bands 
on three occasions (see Table In run 1 we observed 
the star during twelve nights with the 0.91m ESO-Dutch 
telescope at La Silla. The observing schedule was divided 
in four blocks of three nights; every first night the star 
was observed for an average of 5 consecutive hours in a 
U BV Gunn i-exposure sequence while every second and 
third night typically one to three exposures were taken 
in each filter. During both nights of run 2, S 1082 was 
observed for several hours with the lm Jacobus Kapteyn 
Telescope on La Palma. In run 3, also on the Jacobus 
Kapteyn Telescope, we aimed to complete the phase cov- 
erage of the light curve between phases 0-0.3. 

Standard reduction steps of bias subtraction and flat- 
fielding were performed with IRAF[] routines. Aperture 
photometry for all the stars in the field was done with the 
DAOPHOT.PHOT task. Differential light curves for each in- 
dividual dataset were computed with ensemble photom- 
etry (Honeycutt 1992). The variability properties of the 
other stars in the fields are discussed in van den Berg et 
al. (2001) and Stassun et al. (2001, in preparation). We 
refer the readers to these papers for a full description of 
the observations and the photometry reduction. 

2.2. Spectroscopy 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation 



2.2.1. High-resolution spectra 

High-resolution (R«49 000) echelle spectra were taken 
with the Utrecht Echelle Spectrograph on the 4.2m 
William Herschel Telescope on La Palma. S 1082 was ob- 
served on four nights in 1996 and 2000 (see Table || for a 
log of all spectroscopic observations) . 

The 1996-spectra were centred on a blue (4250 A) and 
red (5930 A) wavelength. The 31 lines mm grating was 
used in combination with the 1024x1024 pixels 2 TEK- 
CCD. For a full description of the spectra of run 1 we 
refer to van den Berg et al. (1999). 

The spectra of 2000 were all taken with the same in- 
strumental setup: the 79 lines mm -1 grating was used with 
the 2148x2148 pixels 2 SITel-CCD while the spectra were 
centred on 5584 A. In run 2 the seeing was about 2", while 
light clouds were present during the start of the run. The 
slit width was set to 1". To secure stability, no changes 
were made to the instrumental setup during the night. In 
the first two observations of run 5 the seeing was about 
2"; this deteriorated to 3" with cloudiness during the last 
two observations. As the slit width was kept fixed at 1", 
these spectra are of bad quality. In run 6 the seeing was 
2" during the first two observations but later improved to 
l'/6. The slit width was accordingly changed from 2 to 1'.'2. 
Due to the wider slit these spectra have a lower resolution 
than the spectra of run 1 and 5. All frames were exposed 
for 1200 s, except for those of run 5 that were exposed for 
1800 s to account for the bad seeing conditions. During 
each run we observed radial- velocity standards. Flatfield 
images were made with exposures of a Tungsten lamp. 
Thorium- Argon lamp emission-line spectra were taken for 
wavelength calibration. 

Data-reduction was done in IRAF with CCDRED and 
echelle routines. After correcting the frames for the elec- 
tronic bias and after flat fielding, spectra for each echelle- 
order were extracted with optimal extraction. Towards 
the red, gaps occur in the wavelength coverage. For the 
wavelength calibration, fifth-order polynomials were fit- 
ted in both directions to the positions of the arclines on 
the CCD; the maximum residuals to the fit were ~0.0125 
A corresponding to 0.75 km s -1 . Low-order polynomials 
were fitted to the spectra for continuum normalization. 



2.2.2. Intermediate-resolution spectra 

Two intermediate-resolution (R«3600) spectra were ob- 
tained with the Intermediate Dispersion Spectrograph 
(IDS) mounted on the 2.5m Isaac Newton Telescope on 
La Palma, on February 20.90214, 2000 (UT) and February 
22.85523 (UT). The R1200B grating and the EEV10 CCD 
combination gave a spectral resolution of 1.19 A (FWHM) 
and a useful wavelength range of 3533-4825A. The expo- 
sure times were 900 seconds each, and the signal-to-noise 
level in the extracted spectra ranged from about 65 per 
pixel near the Balmer jump to about 150 near 4750 A. 
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Table 2. Log of the spectroscopic observations. For each 
run we give the run number, the date of the observations, 
the wavelength coverage in A, the number of observations, 
the exposure time in seconds, the spectrograph used and 
(if applicable) the radial- velocity standards. 



Run Dates 


AA 


# 


tcxp 


Inst 


RV 


1 


Feb 28 1996 


3920-4920 


2 


600, 


ues 


HD 136202 










300 










4890-7940 


3 


360 


ues 


HD 136202 


2 


Feb 16 2000 


4380-8650 


17 


1200 


ues 


HD89449 6 


3 


Feb 20 2000 


3535-5035 


1 


900 


ids 




4 


Feb 22 2000 


3535-5035 


1 


900 


ids 




5 


Mar 13 2000 


4380-8650 


4 


1800 


ues 


HD 89449 


6 


Mar 20 2000 


4380-8650 


5 


1200 


ues 


HD 89449 



a F8III-IV, 6 F6IV (Simbad) 
3. Results 

3.1. Eclipse light curve 

The new light curves of S 1082 confirm the findings of 
Goranskij et al. (1992) and fill in the gap in their data 
between photometric phase 0.53 and 0.74 (with phase 
corresponding to primary minimum). In Fig. [I] our data 
are folded on their ephemeris for the primary minimum: 



Mini = 2 444 643.253(5)+ 1^0677978(50) E 



(1) 



Two unequal eclipses with an amplitude of about 0.09 
and 0.05 mag in V are clearly visible with the deeper, pri- 
mary eclipses occuring around phase as expected from 
Eq. [|. Goranskij et al. note that the light curve near sec- 
ond quadrature (phase 0.75) is systematically higher than 
near first quadrature (phase 0.25) by 0.01 to 0.02 mag in 
V. This might also be the case in our data, but the scatter 
between phases 0.1 and 0.3 makes this difficult to see. The 
scatter is probably related to the bad observing conditions 
(nearby moon) of the last run. 

3.2. Spectral line profiles 

The spectra of S 1082 are dominated by a narrow-lined 
component resembling a star of type early-F. In addition, 
the presence of a broad and shallow component in the 
high-resolution spectra is prominent in the Ha line (see 
also van den Berg et al. 1999). When the spectra are 
arranged in order of light curve phase using Eq. |], it is 
clear that the position and strength of the asymmetric 
absorption in the wings vary regularly (Fig. ||); the spec- 
tra that were included in Fig. 8 of van den Berg et al. 
(1999) correspond to photometric phases 0.092, 0.16 and 
0.25. Around first quadrature (phase 0.25) the depth of 
Ha is smallest while the asymmetric absorption in the 
wings is maximally blue-shifted; for phases between 0.5 
and 1 the phase-coverage is sparser, but in the spectra of 
phases 0.88 and 0.90 it can be seen that the asymmetric 
absorption has moved to the red wing. The timing of the 
velocity shift with respect to the eclipses associates the 
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Fig. 1. Light and colour curves for S 1082 folded on the 
ephemeris of Goranskij et al. (1992). Data from the three 
runs are indicated with different symbols: open circles for 
run 1, filled circles for run 2 and triangles for run 3. A typi- 
cal error bar is shown for each run. Magnitude and colours 
are plotted relative to the average V (11.25), U-V (0.45), 
B-V (0.41) and V—I (0.57) values in our measurements. 



broad-lined feature with the brighter star in the eclipsing 
binary (the brighter star must approach us after the pri- 
mary eclipse). Similar variable line profiles are visible in 
the HP, Nal D, Mgl b and Call infrared lines (Fig. |). 

The temperature of the narrow-lined star in S 1082 
can be measured with the low-resolution spectra. We as- 
sume that its light is least contaminated by its companion 
star(s) in the blue and therefore concentrate on the region 
around the Balmer jump. In Fig. ^ synthetic spectra for 
various effective temperatures T e g are compared with the 
IDS spectrum between 3535 A and 4050 A. The spectra 
were computed with model atmospheres for solar metallic- 
ity (Kurucz 1979). The observed spectrum was corrected 
for the reddening towards M 67 (E(B - F)=0.032±0.006, 
Nissen et al. 1987) and for a radial-velocity difference with 
respect to the model spectra. From the relative strength of 
the Call H&K lines, a sensitive temperature indicator in 
this region (e.g. Gray & Garrison 1989), and the contrast 
in the Balmer jump it is clear that the observed spectrum 
is hotter than that of a 6500 K star. 
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Fig. 2. Series of Ha profiles in S 1082. The phase of the 
light curve is indicated for each spectrum. Note the central 
depth of the line and the asymmetric wings. The observa- 
tions during phase 0.54 and 0.58 are noisy due to the bad 
weather conditions of run 5. The continuum separation 
between the spectra is 0.5 unit. 



This part of the spectrum was fitted to model spectra 
with surface gravity logg ranging from 0.5 to 5.0 in steps 
of 0.5, T cS from 6000 to 8000 K in steps of 250 K and a 
fixed projected rotational velocity v TO t smi of 10 km s _1 . 
A free parameter is a wavelength-independent scale factor 
ranging from 0.025 to 1 in increments of 0.025 that is a 
measure of the relative contribution of the narrow-lined 
star. Observed and model spectra are normalised to the 
flux at 4050 A. A straight line is fitted to the difference be- 
tween the observed and each scaled model spectrum. The 
model that produces the smallest residuals to the fit has 
T ef j=7500 K, log <7=4.5 and scale factor 0.85. If the tem- 
perature of the hot star is indeed 7500 K, it is a late A 
rather than an F star; the value of log g is close to that of 
a main-sequence star for which log g w 4.25 (Gray 1992). 
This shows that the narrow-lined component dominates 
the spectrum in the blue and it implies that one of the in- 
dividual components of S 1082 remains a blue straggler. A 
more accurate decomposition of S 1082 is given in Sect. |]. 
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Fig. 3. H/3, Mglb and Call IRT lines in S 1082 for a se- 
lected set of observations. Light curve phases are indi- 
cated. The continuum separation between the spectra is 
0.5 unit. 



3.3. Radial-velocity curves 

Radial velocities were measured via cross correlation of the 
high-resolution spectra with appropriate template spec- 
tra (Table ||). To determine the radial velocities of the 
narrow-lined star in S 1082 we adopted as templates spec- 
tra of F-type radial- velocity standards observed during the 
same run (Table §). Velocities were derived for each order 
separately. Only orders without strong telluric lines were 
selected, and included 4890-6820 A for the 1996-spectra 
and 4435-6820 A for the 2000-spectra. The radial veloci- 
ties listed in the third column of Table |^ are the averages 
of the values from the individual orders; the errors repre- 
sent the spread around the average. Note that systematic 
errors can still be included, e.g. due to the wavelength 
calibration. We exp ect the latter not to exceed 0.75 km 
s _1 (see Sect. 2.2. l|) . Our measurements confirm that the 
narrow-lined component in S 1082 shows radial-velocity 
variations of only a few km s _1 . 

The radial- velocity curves of the stars in the eclipsing 
binary were measured after the contribution of the hot, 
narrow-lined star to the total spectrum was removed. To 
that end, all spectra were first brought to the rest frame of 
the hot star. Then, a hot-star template spectrum was con- 
structed by taking the median of the spectra obtained at 
phases 0.95, 0.97, 0.0035, 0.019 and 0.048; these particular 
phases were chosen in order to let the spectral line pro- 
files of the template be as symmetric as possible. For the 
1996-spectra, taken with a different instrumental setup, 
the spectrum obtained at phase 0.092 was chosen as a 
template. The relative contribution of the hot-star tem- 
plate to the remaining spectra was derived with the same 
fitting procedure as described in Sect. ^7J: the template 
multiplied with scaling factors ranging from 0.025 to 1 
was subtracted from each individual spectrum, and the 
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wavelength (A) 

Fig. 4. Spectrum of S 1082 compared with synthetic spec- 
tra of various effective temperatures; surface gravity and 
projected rotational velocity are fixed at logg=4.5 and 
D ro tsinz=10 km s~ 1 . The spectra are normalised to the 
flux at 4050 A; each spectrum is offset with 1 unit from 
the previous spectrum. The relative strength of the Call 
K line and the Call H+He blend (marked with 'K' and 
'H'), and the contrast in the Balmer jump (flux near 4050 
A relative to flux near 3600 A) are indicators of tempera- 
ture. 

scaling that produces the smallest residuals around a fit 
to a straight line was chosen as the appropriate scaling for 
that particular spectrum. The scaled template was then 
subtracted from the total spectrum to obtain the spectrum 
of the binary at various phases. 

Next, these residual spectra were correlated against a 
synthetic spectrum of T e ff=5250 K and log g=3.5. For this 
model spectrum we choose z; ro t sinj=50 km s _1 to roughly 
match the apparently broader lines in the spectrum of the 
binary. This procedure was repeated for every order be- 
tween 4380 and 6435 A, excluding the order that contains 
the broad H/3 line and the region between 5690 and 6090 



Table 3. Radial velocities of the three stars in S 1082, the 
components of the eclipsing binary Aa and Ab, and the 
outer companion B. From left to right: heliocentric Julian 
date (-2 450 000) at the midpoint of observation; photo- 
metric phase computed with the ephemeris of Goranskij 
et al. (1992); heliocentric radial velocity for star B; radial 
velocity of the primary Aa and secondary Ab in the eclips- 
ing binary with respect to star B in km s _1 . The spectra of 
phases 0.5378 and 0.5609 were taken under bad observing 
conditions. 



riJJJ- 


phot. 


frad,B 


^rad,Aa 


«rad,Ab 


I 4oU UUU 


phase 


km s — 


km s~ 1 


km s — 


1 A O CI Al 


nooi o 
.09213 


oo n_i_n 1 
33.U±U.1 






i A o c 7nn 
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oo o_i_n i 
33.2±U.l 


n^-L a 
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142.6837 


.2548 


oo n I n o 

32. 9±U. 2 
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i Km oooo 
to91.3o22 


n'ynn 
.9/09 


OO 7_LA O 
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t591.4t lU 


nno a k i 
.003451 


oa o_i_n o 






i Km aoo a 
t591.4334 


m oon 
.01850 


oo 7xn o 
33. (±0.3 






i Km Ad An 
lt>91.464U 


.04/50 


oo a _i_n o 






i Km a one 
1591. 4505 


neon a 
.U6294 


oo o_i_n o 
32. 8±U. 3 






1591.4967 


.07812 


32.7±0.3 


-73±4 


122±7 


1591.5452 


.1235 


33.2±0.4 


-96±3 


148±4 


1591.5616 


.1389 


33.3±0.4 


-95±2 


141±4 


1591.5779 


.1541 


33.3±0.4 


-100±2 


153±5 


1591.5993 


.1742 


33.1±0.5 


-104±2 


154±5 


1591.6150 


.1890 


33.5±0.4 


-108±2 


150±12 


1591.6310 


.2039 


33.9±0.5 


-107±3 


166±5 


1591.6474 


.2192 


34.2±0.5 


-109±2 


170±8 


1591.6682 


.2387 


33.8±0.5 


-113±2 


160±7 


1591.6840 


.2535 


33.9±0.5 


-113±3 


169±14 


1591.7010 


.2695 


33.9±0.5 


-115±2 


180±5 


1617.4072 


.3435 


31.8±0.4 


-104±2 


163±5 


1617.4308 


.3656 


32.7±0.5 


-96±2 


153±5 


1617.6147 


.5378 


36. Oil. 2 


28±5 




1617.6394 


.5609 


35.2±0.8 






1624.3903 


.8832 


35.4±0.4 


91±3 


-130±5 


1624.4092 


.9009 


34.6±0.4 


85±4 


-132±7 


1624.5766 


.05765 


32.4±0.3 






1624.5949 


.07484 


32.7±0.3 


-80±4 


121±6 


1624.6131 


.09188 


32.6±0.4 


-90±4 


128±5 



for which the cross correlation functions are very noisy 
with no clear peaks. A cross correlation peak is some- 
times visible at w ra d=0 km s _1 and represents features of 
the hot star that were not corrected for (or introduced 
into the spectrum) by subtraction of the template. The 
cross correlation functions of most spectra clearly show 
two peaks with a variable separation (see Fig. |5|). One peak 
is smaller and broader than the other, and is redshifted 
with respect to the template for phases smaller than 0.5; 
therefore, this peak is identified with the secondary star 
in the eclipsing binary. This is the first spectroscopic ev- 
idence for the third star in S 1082. The measurements of 
the velocity of the primary near phase 0.1 are probably 
distorted due to its eclipse by the secondary star. In the 
following we will refer to the components of the eclipsing 
binary as Aa and Ab, and to the outer companion as B. 
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Fig. 5. Cross correlation function resulting from the resid- 
ual spectrum taken at photometric phase 0.2535 after the 
contribution of the hot star is removed. The left peak cor- 
responds to the primary star Aa in the eclipsing binary. 
The width of peak corresponding to the secondary Ab is 
broader which indicates that its spectral lines are more 
broadened by rotation. For this particular cross correla- 
tion the spectrum between 5160 and 5240 A was used. 
The radial velocity is given with respect to the hot star B 
in S1082. 

We use the spectrum observed near first quadrature 
(wavelength region 5160 to 5240 A containing the Mgl b 
triplet) to derive an estimate for the projected rotational 
velocities of both components from the widths of their 
cross correlation peaks. To that end, we cross correlate a 
synthetic spectrum of a certain temperature (logg = 4.0, 
u ro t sini=30 km s _1 ) with synthetic templates of the same 
temperature and surface gravity with w ro t sin i ranging 
from 10 to 100 km s -1 , and with the object spectrum. 
The peak of the cross correlation functions are fitted with 
gaussian profiles. This is repeated for synthetic spectra of 
temperatures 5000 to 6750 K in steps of 250 K. With the 
results we construct calibration curves for each tempera- 
ture that give the width of the cross correlation peak as 
a function of the projected rotational velocity of the syn- 
thetic template. The widths of the cross correlation peaks 
from the object spectrum then give a rough estimate for 
the w ro t sin i of both stars. For star Aa we obtain a v m t sin i 
of 51 km s- 1 (5750 K) to 58 km s" 1 (6750 K), for star Ab 
a v mt sin i of 83 km s" 1 (5500, 5750 K) to 87 km s" 1 (6750 
K). 

4. Parameter estimation for the eclipsing binary 

We used the recently written Eclipsing Light Curve (ELC) 
code of Orosz & Hauschildt (2000) to fit the basic binary 
observables of S 1082, namely the light curves in the four 
filters, the two velocity curves and the observed projected 
rotational velocities of the two stars Aa and Ab. The blue 
straggler B in this case is modelled as "third light" . The 
addition of a constant (in phase) light will reduce the ob- 



served relative amplitudes of the light curves. The ELC 
code uses model atmosphere specific intensities, so the flux 
from a third star can be self-consistently added to several 
different bandpasses (i.e. a blue third light star dilutes the 
5-band light curve of a red binary more than it dilutes the 
/-band light curve). We mainly used the specific intensi- 
ties derived from the NextGen models (Hauschildt et 
al. 1999a; 1999b). A few models were also computed us- 
ing specific intensities derived from Kurucz (1979) models. 
Only the photometric data of run 1 were used for the light 
curve fitting as they form the largest set of homogeneous 
measurements. 

We start by estimating the component masses of the 
eclipsing binary. Unfortunately, we only have radial veloc- 
ities near one quadrature phase (poor weather prevented 
us from observing the opposite quadrature), so the veloc- 
ity amplitudes of the two curves and the binary systemic 
velocity will not be as constrained as we would like them 
to be. We fit sinusoids to each of the velocity curves, fix- 
ing the period at the photometric period. For star Aa we 
find K Aa = 120.8 ± 6.4 km s~\ 7 rel)Aa = 1.7 ± 5.8 km 
s _1 , 0o,Aa(spect) = 0.747 ± 0.007, where 7 ro i is the sys- 
temic velocity relative to the mean cluster velocity and 
where </>o(spect) refers to the phase of maximum veloc- 
ity. For star Ab we find K A b = 190.1 ± 12.5 km s -1 , 
7rel,Ab = -2.6il2.5kms" 1 , (£o,Ab (spect) = 0.242±0.010. 
The errors on the individual velocities were scaled to yield 
Xv = 1 for each curve, and the error estimates on the fit- 
ted parameters were derived using the scaled uncertain- 
ties. The phasing of the curves are consistent with ex- 
pectations, where star Aa has its maximum velocity one 
fourth of an orbital cycle after the deeper photometric 
eclipse. Since systemic velocities are consistent with be- 
ing zero, we will assume the binary is a cluster member 
and hence has 7 rc i = km s . In this case we find for 
star Aa i^Aa = 119.1 ± 3.1 km s _1 and 0o,Aa(spect) = 
0.748 ± 0.007, and for star Ab we find K A b = 190.1 ± 11.6 
km s" 1 and o ,Ab(spect) = 0.243 ± 0.010. Taking the 
sinusoid fits at face value, we can immediately compute 
the mass ratio of the binary. We find Q = Ma^/Mah = 
KAa/KAb — 0.63 ± 0.04. The minimum masses of the two 
component stars are then Ma s sin 3 i — 2.01 ±0.38 Mq and 
M Ab sin 3 i = 1.26 ± 0.27 M (see Table || for a summary 
of the observed parameters of the eclipsing binary) . 

There are a total of twelve free parameters in the 
model: the filling factors (by radius) of the two stars /A a , 
/Ab) the mean temperatures of the two stars TA a , ?Ab> 
the "spin factors" of the two stars f^Aa, ^Ab, where Q is 
the ratio of the rotational angular velocity to the orbital 
angular velocity, the inclination i, the mass ratio Q, the 
orbital separation a, the temperature of the third light 
star Tb, the surface gravity of the third light star loggs, 
and the third light scaling factor SAb- We assume the 
orbit is circular, and that the third light star B does not 
vary. The gravity darkening exponent of star Aa was fixed 
at 0.25, the standard value for a star with a radiative en- 
velope, while the gravity darkening exponent for star Ab 
was set at the standard convective value of 0.08. The ELC 
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Fig. 6. Fit to the light curves in the U, B, V and Gunn i band. The small panels show the residuals to the fit. Only 
the data from the first observation run are used. 



code uses Wilson's (1990) detailed reflection scheme, and 
for this problem the albedo of star Aa was taken to be 
1 and the albedo of star Ab was taken to be 0.5. Three 
iterations of the reflection scheme were needed to achieve 
convergence. Both stars in the binary are assumed to be 
free of spots. A variation of the "grid search" routine out- 
lined in Bevington (1969) was used to optimise the fits. 
In practice the fitting procedure involved a great deal of 
interaction where some parameters were temporarily fixed 
at certain values. Several two-dimensional grids in param- 
eter space were defined (for example a grid of points in the 
/Aa, /Ab plane). For each point, we fixed those parameters 
at the values defined by the grid location and optimised 
the other parameters, creating contours of \ 2 values. New 
parameter sets were optimised using the set of parameters 



for a nearby point that gave the best previous fit. After the 
lowest x 2 value in a grid was found, then a new grid using 
other parameters was computed using the best solution 
as a starting point. The fitting took several weeks of CPU 
time on an Alpha XP 1000, and sampled a wide range 
of parameter space. We are reasonably confident that our 
results are at or very near the global \ 2 minimum. 

We found a relatively large number of solutions with 
similar \ 2 values. Fig. ^ shows a typical fit. The light 
curves in the 4 bands (Fig. ^) are fitted reasonably well, 
although there are still some systematic deviations, espe- 
cially near the eclipse phases. The radial velocities are also 
fitted reasonably well (Fig. 0), but again there are some 
small systematic deviations (the velocity curves seem to 
be systematically flatter than the sinusoid fits near phase 
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0.25). In all solutions the binary is detached, i.e. both stars 
are well within their respective Roche lobes. We applied 
two additional constraints in order to narrow down the 
range of parameters. The first constraint is that the to- 
tal V magnitude and B — V colour of the model should 
match the observed V magnitude of S 1082. In this case 
the apparent V magnitude of the model is easy to com- 
pute. We used the synthetic photometry computed from 
the NextGen models^] to compute expected absolute V 
magnitude of each star Aa and Ab in the eclipsing binary 
from its temperature, radius, and surface gravity. The V 
magnitude of the blue straggler B then follows from the 
fitted luminosity scaling. The second constraint is that 
the implied mass of the blue straggler is roughly consis- 
tent with its place in the colour-magnitude diagram. That 
is, the radius of the blue straggler B can be computed 
from the distance, the V magnitude, and the temperature. 
Since the gravity of the blue straggler B is specified in the 
models, its mass can then be computed. The blue colour of 
S 1082 requires a relatively hot third light star (« 7500 K), 
and its surface gravity must be near logg = 4.25 in order 
for the mass to be near ps 1.7 Mq. The derived astrophys- 
ical parameters for the adopted model are summarised in 
Table §. The errors on the parameters were estimated from 
the x 2 values generated in the various grid searches. These 
error estimates may be too small, given the complicated 
nature of the model. The errors on the masses were taken 
to be on the order of 15%, as judged from the quality of 
the sinusoid fits. The light curves of the three components 
are shown in Fig. ||, and a cartoon of the binary at three 
phases in Fig. [| 

We have assumed that the stars are not spotted and 
that the eclipsing binary has a circular orbit. A violation 
of either one of these assumptions could alter the light 
curves to produce the small systematic deviations seen in 
the residuals. Bright or dark spots could either add or 
remove light at certain phases, complicating the analysis. 
Given the rather large number of free parameters we have 
now we did not consider adding spots at this time. If this 
binary is part of a triple, then the orbit could be eccentric 
(see Sect. ||). A slight eccentricity (e « 0.05 say) could 
cause the maxima to be asymmetric and the minima to 
be shifted slightly in phase. Our current velocity curves 
do not have enough phase coverage to place meaningful 
constraints on the eccentricity, so any firm conclusions on 
the eccentricity will have to await the arrival of additional 
data. 

Fig. [To| shows the decomposition of S 1082 into its com- 
ponents in a colour-magnitude diagram of M 67. Star Ab is 
located near the 4 Gyr isochrone that provides the best fit 
to the observed main-sequence, subgiant and giant stars 
(Pols et al. 1998). Hence the position of star Ab is con- 
sistent with the expected value based on its mass. The 
positions of star Aa and the blue straggler B are a bit 
more uncertain since one can "trade off" flux between the 
two stars (i.e. the blue straggler B can be made brighter 

2 [ftp: //calvin.physast .uga. edu/pub/NextGen/Colors/ 



Table 4. Observed binary parameters. Aa and Ab refer 
to the stars in the eclipsing binary. V and B — V are taken 
from Montgomery et al. (1993). 



parameter 


value 


-Pphot 


1.0677971(7) 


To 


2 444 643.250(2) 


V 


11.251 


B-V 


0.415 


Ka& (km s -1 ) 


119.1(3.1) 


K A b (km s _1 ) 


190.1(11.6) 


Q 


0.63(4) 


7Aa (km s 1 ) 


0.0 [fixed] 


7Ab (km s -1 ) 


0.0 [fixed] 


M A asin 3 i (Mq) 


2.01(38) 


A/ Ab sin s i (Mq) 


1.26(27) 


u ro t,Aa sini (km s _1 ) 


56(5) 


«rot,Ab sini (km s _1 ) 


83(5) 



Table 5. Fitted binary parameters. Aa and Ab refer to 
the stars in the eclipsing binary, B to the third, outer star. 



parameter 


Aa 


Ab 


B 


/ 


0.520(10) 


0.700(10) 




R (Re) 


2.07(7) 


2.17(3) 


« 2.5 


R/ -RRochc 


0.66(7) 


0.86(7) 






0.49(5) 


0.91(5) 




T(K) 


6480(25) 


5450(40) 


7500(50) 


log g (cgs) 


4.21(2) 


4.0(2) 


4.25(5) 


M (Mq) 


2.70(38) 


1.70(27) 


« 1.7 


V 


12.33(11) 


13.10(11) 


12.24(11) 


B — V 


0.51(2) 


0.82(2) 


0.33(1) 


ia (deg) 


64.0(1.0) 






a a (Re) 


7.2(4) 






V, B - V total 


11.30(11), 


0.48(2) 





at the expense of star Aa in the binary). We note that 
the star Ab is always located (within the errors) on the 4 
Gyr isochrone. The error bars shown on the V magnitudes 
reflect the uncertainties for our adopted model which pro- 
duces an overall V magnitude of all three stars close to the 
observed value. In any event, the hotter star Aa is sublu- 
minous by at least 2 mag in V (compare with the isochrone 
for a 2.2 Mq star). The position of the blue straggler lies 
to the red of the extension of the main-sequence as defined 
by the 1 Myr isochrone. 

In order to obtain observed times of primary mini- 
mum we use the model light curve as a template to fit the 
data near the primary minima observed during run 1 and 
run 2. Seven additional times of minimum were obtained 
from fitting the data in Table 3 of Goranskij et al. and 
from the data points of Simoda (1991) (see first column of 
Table |6|) . A straight line is fitted to the observed times of 
minimum to find a new period and T (see Table |). The 
period thus derived is compatible with the period listed by 
Goranskij et al. (Eq. [l]). We use the new ephemeris to com- 
pute observed minus computed (O-C) times of primary 
eclipse and the corresponding cycle number with respect 
to T (second and third column of Table |]). The peak-to- 
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Fig. 7. Radial-velocity curves of star Aa of the eclipsing 
binary in S 1082 (filled circles) and star Ab (filled trian- 
gles). The lines indicate the best- fitting velocity curves as 
computed by ELC. 
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Fig. 8. Light curves of the three components in S 1082 
(dashed-dotted line indicates the primary Aa in the eclips- 
ing binary). 



peak amplitude AO-C is ~39 minutes. If these variations 
are real and caused by the motion of the eclipsing binary 
around a third body it would correspond to a minimum 
semi-major axis in the outer (o) orbit of the binary (b) 
a Q ,b sini =l/2 c AO-C=2.3 AU. Assuming that the blue 
straggler (B) has a mass compatible with its position in 
the colour- magnitude diagram (about 1.7 M Q ) this corre- 
sponds to a minimum value of the semi-major axis of the 
total system a sini =a 0j b sini (1+Ma/A/b) ~8.4 AU; 
combined with the total mass of the system of 6.1 Mq and 
Kepler's third law this gives a minimum period of 10 years. 
This is not compatible with the period of ^3 years found 
by Milone (1991). We conclude that not all the variation 
in O-C is due to light-travel time in the outer orbit. There 
is no evidence for periodicity in the O-C times, although 
the time baseline is somewhat short and the coverage is 
somewhat spotty. 




Fig. 9. Cartoon of the configuration at three orbital 
phases of the inner binary according to the model listed 
in Table || Note that the Roche-lobes of the two stars Aa 
and Ab slightly overlap due to the fact that neither star 
corotates with the orbit. 

Table 6. Times of primary minimum. From left to right: 
observed time of minimum; O-C difference (in days) be- 
tween the observed and computed time of eclipse; cycle 
number of the eclipse with respect to the To of Table |[ 
The errors in the observed times of minimum and the O-C 
values are 0.005 days. 



HJD 


O-C 


cycle 


source 


-2 440 000 








3191.036 


-0.010 


-1360 


Simoda (1991)° 


4643.250 


0.000 





Goranskij et al. (1992) 


5325.586 


0.013 


639 


idem 


6773.492 


-0.013 


1995 


idem 


7861.605 


0.014 


3014 


idem 


7920.333 


0.014 


3069 


idem 


7944.869 


-0.010 


3092 


idem 6 


11218.744 


-0.001 


6158 


run 1 


11539.078 


-0.006 


6458 


run 2 



"Simoda observed two consecutive primary eclipses. These 
data were combined to measure one time of minimum. 
6 Based on a measurement of a secondary eclipse which we con- 
vert to a primary eclipse by adding half a period. We note that 
if the orbit of the eclipsing binary is eccentric, the two eclipses 
need not be separated by half a period. 

5. Discussion 

Our new observations of S 1082 confirm the eclipses re- 
ported by Goranskij et al. (1992), the small radial- velocity 
variations of the narrow lines in the spectrum (Mathieu et 



10 Maureen van den Berg et al.: The blue straggler S 1082: a triple system in the old open cluster M67 





B-V 

Fig. 10. Colour-magnitude diagram of M67 that shows the decomposition of S 1082 into a blue straggler (B) and 
stars Aa and Ab in the eclipsing binary (A); a dashed line connects their positions. The observed location of S 1082 is 
indicated with a box; the height of the square equals the depth of the primary eclipse. The dotted lines arc evolutionary 
tracks for 1.6, 1.8, 2.0 and 2.2 M Q stars corrected for the distance modulus and reddening of M67 (.2=0.02, Pols et 
al. 1998). The 4 Gyr-isochrone is included to indicate the expected positions of cluster members; the 1 Myr isochronc 
is included to give an estimate for the location of the ZAMS (Pols et al. 1998). B and V magnitudes of M67 stars are 
from Montgomery et al. (1993). Only stars with a proper-motion membership probability > 0.8 (Girard et al. 1989) 
are plotted. 



al. 1986) and the time-variation of a broad-lined compo- 
nent seen in high-resolution spectra (van den Berg et al. 
1999, Shetrone & Sandquist 2000). In addition, we now 
have found that the radial-velocity shifts of the broad- 
lined component and of a third component in the spectra 
vary on the photometric period. This clearly demonstrates 
that the narrow-lined star, which is a blue straggler on it 
own, is not part of the eclipsing binary - this solves the 
seeming contradiction in the properties of this system. The 



broad variable component of the Ha line discussed by van 
den Berg et al. (1999) is the Stark-broadened Ha line of 
the hot component Aa of the inner binary. The ultravio- 
let flux measured by Landsman et al. (1998), higher than 
expected for a star at the B — V colour of S 1082, is only 
partly explained as a consequence of the bluer colour, i.e. 
higher temperature of the third star. 

Comparison with Figs. 2 and 3 of Dempsey et al. 
(1993) shows that the X-ray luminosity of S 1082 (Belloni 
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et al. 1998) is typical for a subgiant in an RS CVn system. 
The X-ray emission could therefore be caused by mag- 
netic activity in the rapidly rotating subgiant Ab. Thus, 
our model can explain the X-rays of S 1082. 

Our solution to the photometry is symmetric (see 
Fig. ^|) and does not explain the asymmetry between 
phases 0.25 and 0.75; nor the variability in the form of 
the secondary eclipse between our first two observations 
runs. This indicates that variable spots are present and 
accordingly that the best-fit parameters of the system are 
subject to some additional uncertainty. However, we do 
not think our main conclusions are affected. 

We think that the eclipsing binary forms a bound triple 
with the third star, for two reasons. First, the radial- 
velocity measurements of the third star B - the narrow- 
lined system in S 1082 - indicate that it is in a ~1000day 
orbit around a companion. The systemic velocities of both 
the 1000 day and the 1.07 day binaries are compatible with 
the radial velocity of the cluster. A chance alignment of 
two such bright cluster members is unlikely. Second, if we 
relax the constraint that the binary is at the cluster dis- 
tance in solving the photometry, we find that the best solu- 
tions tend to be an Algol binary at higher mass (> 4 Mq) 
and larger distance (about twice the cluster distance). A 
high-mass binary at such a large distance from the galactic 
disk is unlikely. 

Assuming then that the binary forms a hierarchical 
triple with the third star, in M67, we note that the 1-er 
lower bound to the mass of the binary is at about three 
times the turnoff mass of the cluster. The formation of 
the binary must thus have involved at least three stars. 
The third star is a blue straggler on its own account, and 
thus according to most current models its formation in- 
volves two stars. We have to conclude that the formation 
of S 1082 required the interaction of no less than five stars! 

This interaction may have started with a binary-binary 
encounter, in which two stars collided directly and merged; 
one of the remaining two stars ended in a close orbit 
around the merger, the fourth star in a wide orbit around 
the inner binary. The fourth star must be a blue straggler, 
which either was already present in the original encounter 
or was later exchanged into the system. That triples are 
formed easily by binary-binary encounters, and that they 
live long enough to undergo subsequent exchange encoun- 
ters with a binary, is shown for example by the compu- 
tations of Aarseth & Mardling (2001). It is less obvious 
that mergers are common in such encounters, as most 
binary-binary encounters are between relatively wide sys- 
tems. Merger products tend to be subluminous (Sills et al. 
2001), in accordance with the properties of the more mas- 
sive star Aa in the inner binary. The main problem with 
this scenario for the formation of the S 1082 system is that 
a merger product is only subluminous for a very limited 
period; comparable to its thermal time scale. This reduces 
the available time within which the outer star of the first 
encounter is exchanged with a binary or blue straggler in a 
subsequent encounter, and indicates that the blue strag- 
gler was present in the initial binary-binary encounter. 



Alternatively, an initial encounter between two triple sys- 
tems can have led immediately to the currently observed 
configuration. In any case, the probability of catching the 
merger product while it is strongly subluminous is discon- 
certingly small. 

Aarseth & Mardling (2001) note that the inclination 
of the outer orbit in a hierarchical triple with respect to 
the inner orbit can induce a large eccentricity in the inner 
orbit; subsequently, tidal forces in the inner orbit cause 
it to shrink. Thus, the inner orbit may be eccentric, and 
smaller than in the past. Our solution for the inner bi- 
nary implies that the cool star has a radius which is a 
sizeable fraction of the binary separation R^h/a-A ~ 0.3. 
According to Eq. 2 of Verbunt & Phinney (1995) the cir- 
cularisation time scale of the inner binary A is ~ 10 3 year; 
the eccentricity of A and asynchronicity of Ab thus will 
be determined by the competing effects of the perturba- 
tion by the outer star B and the tidal forces in the inner 
binary. Only if the latter are substantially reduced with 
respect to the Zahn (1977) formulation - as suggested by 
e.g. Goodman & Oh (1997) - do we expect measurable 
eccentricity and asynchronicity. 

Whereas our new observations have allowed us to re- 
solve the apparent contradictions in the earlier data and to 
determine the system parameters, we conclude that these 
parameters - in particular those of the primary Aa of the 
inner binary - are hard to understand in terms of standard 
stellar and binary evolution, even when stellar encounters 
and mergers are taken into account. This situation is re- 
markably similar to that in the study of two other mem- 
bers of M67, located below the subgiant branch in the 
colour-magnitude diagram, both of which turn out to be 
binaries (Mathieu et al. 2001, in preparation); and indeed 
to our lack of understanding of the blue straggler popula- 
tion in the cluster. Many of these stars may be mergers, 
which raises the question whether the object that results 
when two stars merge may follow an evolutionary track 
which is very different from that of an ordinary star of the 
same mass; and whether they can do so for a period of 
time which significantly exceeds the thermal time scale. 

Further studies of the remarkable triple system S 1082 
should provide a better-sampled radial- velocity curve, nec- 
essary to determine whether the inner binary is eccentric, 
and to measure its systemic velocity more accurately to 
establish cluster membership. Better-quality light curves 
may also serve to improve the accuracy of the inner or- 
bit, e.g. if the inner orbit is eccentric the separation of the 
eclipses is different of 0.5. A long-term sampling of the 
light curves will be required to improve the interpretation 
of the O-C in the timings of the primary eclipse. 
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Fig. 11. O-C times of the primary eclipse in days versus 
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